Is turbulent mixing a self-convolution process ? 
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Experimental results for the evolution of the probability distribution function (PDF) of a scalar 
mixed by a turbulence flow in a channel are presented. The sequence of PDF from an initial skewed 
distribution to a sharp Gaussian is found to be non universal. The route toward homogeneization 
depends on the ratio between the cross sections of the dye injector and the channel. In link with 
this observation, advantages, shortcomings and applicability of models for the PDF evolution based 
on a self-convolution mechanisms are discussed. 
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Predicting the temporal evolution of the probability 
distribution function (PDF) of a tracer is of wide inter- 
est, in particular for the dynamics of density stratified 
fluids or for reactive flows [1| . The scope of existing phe- 
nomenological approaches is to provide simple, efficient 
models to describe the evolution of the PDF, from an ar- 
bitrary initial condition to a sharp peak centred around 
the mean, representing the late stage of mixing. 

One of the simplest and widely used models is the lin- 
ear mean square estimate (LMSE) also called interaction 
by exchange with the mean (IEM) [2J. The main limita- 
tion of this model is that the PDF keeps its shape while 
it contracts around the mean by mixing. Many studies 
have been devoted to the improvement of this approach 

sina. 

There exists an alternative class of models, based on a 
self-convolution process. The first occurrence of such a 
model is the coalescence dispersion (CD) mechanism of 
[H , originally intended to describe drop interactions in a 
two liquid system: two drops having independent scalar 
values (Jo and o~ coalesce into one drop having the aver- 
age scalar value a = (oo+^o) which immediately splits 
into two drops with the same scalar value a. The proba- 
bility p to measure a given scalar value after coalescence 
is the convolution of the same probability function po 
before coalescence : p (a) — 2 J po (2a — <tq) po (ctq) da$. 
CD models with self-convolution have been widely stud- 
ied and extended to turbulent flows [9( although some of 
the original underlying hypothesis become questionable 
(lQ| . Following the ideas of the CD approach, a model 
of aggregation of scalar streaks has been presented in 
flUl. A different physical mechanism has been proposed 
in [12j to justify a pure self-convolution process. This 
model describes the temporal evolution for the coarse- 
grained PDF at a given scale I . The key idea is that ran- 
dom straining transfers scalar fluctuations toward scales 
smaller than I , so that the measured fluctuations are aver- 
aged over scale I. Such coarse-grained averaging is some- 
what equivalent to the smoothing effect of diffusivity. 

The question of the PDF evolution of a passive tracer 
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FIG. 1: Experimental set up 



in a turbulent flow has been addressed ex per imentally in 
different configurations (see for instance [13, fl4|). For a 
flow without a mean scalar gradient, convergence from a 
very skewed PDF toward a gaussian shape has been qual- 
itatively observed, although some departure from gaus- 
sianity was noticed to persist even at late stages of mix- 
ing. 

Surprisingly, there seem to be only few comparisons 
between phenomenological models and the laboratory 
experiments presented above. We consider here a 
fluorescent dye introduced in a channel, where it is 
advected by a turbulent flow. PDF are measured at 
different locations along the channel. The PDF becomes 
sharper and sharper around the mean with increasing 
downstream distance from the injector. The sequence 
of PDF, from an initial skewed distribution to a sharp 
gaussian, is found to be non universal, depending on 
the injection system. Our aim is here to discuss the 
capability of models based on a self-convolution process 
to describe those experiments. 

a. Models We consider an idealised case : the scalar 
field and the turbulent flow are supposed to be statisti- 
cally homogeneous in space. We call pi (a, t) the PDF of 
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a scalar value a at time t, measured with a probe of size 
I. The scalar field measured by this probe is averaged at 
scale I by a low pass filter. The self-convolution model 
relies on two hypothesis : 1) it ignores fluctuations of 
the straining histories of scalar sheets 2) it assumes in- 
dependence of scalar values at two points separated by a 
distance /. 

Hypothesis 1) implies that if the width of a scalar 
sheet, equal to I at time t = 0, is divided by 2 at 
time t, then the width of the adjacent sheets are also 
divided by a factor 2. This time interval t is related to 

the strain rate s by the relation exp ^ J Q * s(t')dt'^j = 2. 

The probability to measure a given scalar value at scale 
I at time t is then pi(a,t) — 2 J pi(o\,2o — ai,0)dai, 
where pi(o~\, 02, 0) is the joint probability to measure 
the values o\ and 02 at two points separated by the 
distance I. The hypothesis 2) of independency gives 
pi(a,t) = 2f p l (<7 1 ,0)pi{2a - a 1 ,0)d<7i. 

Taking the limit of continuous times for this process, a 
dynamical equation for the Laplace transform of the PDF 
pi{n) — J +O ° pi(p)e~ Ka do has been provided in (l2T |: 



d t pi = s{t) [pi lnp - nd K pi] 



(1) 



Note that this equation is the continuous form of the 
discrete map used in The solution of ([T]) is 



p(n, t) 



/(f) = exp / s (f) dt' (2) 



From this expression, one can easily deduce (see 
[I2I ]) the temporal evolution of the cumulant c n (t) = 
{-d K ) n \n( Pl {K,t))\ K=0 . c n (t)=c n (0)/f n - 1 . Cumulants 
are equal to the centred moments for 1 < n < 3, and 
related to them for higher n. The application to the sec- 
ond moment yields dc2/dt = —s(t)c2- The cascade rate 
of variance is therefore proportional to the strain rate s. 

The aggregation model provided in leads to a dif- 
ferent dynamical equation: 



d t p = s(t) f(t) 



K,d K p 



(3) 



the solution approaches a sequence of gamma-PDF at 
large time: >y (*/{*)) = J^-ig^e-f 

b. Experimental set up We drive a steady flow in a 
water duct of square section l s xl s = 5 x 5cm 2 , at a steady 
volume rate Q, corresponding to a bulk velocity U = 
Ql (I 2 ). The duct is 150 cm long. Straight circular buffer 
sections, 50 cm long and 5 cm in diameter, are fitted 
at both ends to reduce perturbations from end effects. 
Furthermore, the upstream buffer section is filled by a 
honey-comb. The fluorescent tracer (Rhodamin 6G) is 
introduced at 50 cm downstream from the duct inlet, 
and observed until an abscissa x=80 cm downstream, 
in a region of well established duct flow. Most of the 



experiments have been performed for £7=33 cm/s, so the 
corresponding Reynolds number Re = l s U/v is equal to 
1.6 10 4 . This is about 8 times the Reynolds number of 
transition, so that turbulence is well established. Test 
experiments performed at lower and higher Re (still of 
order 10 ) yield very similar results. 

The mean kinetic energy dissipation rate is estimated 
using standard measurement on turbulent channel flows 
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= (0.316C/ 3 )/(2Z s i?e 1 / 4 ) - 8.10~ 3 m 2 .s 
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The corresponding dissipative Kolmogorov scale is r\ = 
(^ 3 /e) 1 ^ 4 ~ 0.1 mm. To estimate the Taylor microscale 
A = u rms j < (d x u) 2 > 1 / 2 , we consider that turbulent 
fluctuations represent 3.5% of the mean velocity U (fol- 
lowing [l?]]) and use the isotropic relation ^(d x u) 2 ^ = 

e/(15^). This yields A — 0.5 mm, and the corresponding 
Reynolds number is Re\ = u rms \/v ~ 5. This indicates 
the absence of a multiscale energy cascade. 

The Schmidt number of the tracer is Sc = vjD ~ 
1000, corresponding to a very low Batchelor scale lb = 
rjSc^ 1 / 2 ~ 4 pm. Two kinds of dye injectors have been 
used. The first one is a tapped vertical tube with eight 
equidistant small holes, approximating a vertical line 
source with uniform flux. The second is a single tube 
of diameter 2mm, curved in the downstream direction 
with dye flow rate adjusted to minimise shear with the 
background flow. It approximates a point source at the 
centre of the duct. 

Using the Laser Induced Fluorescence technique, the 
scalar field is measured in a vertical plane perpendicu- 
lar to the channel cross section (see figure [T|). The PDF 
of concentration is measured, after proper calibration, 
as a histogram of images from a CCD camera (SMD), 
1024x1024 pixels with 12 bit grey levels. Histograms are 
averaged over series of 300 images, made at a given dis- 
tance x downstream from the injector. We limitate the 
statistics to a central band, 1/4 of the duct width, where 
the mean concentration can be considered as uniform, 
and shearing by the mean flow negligible. Then the vari- 
ation of the PDF with downstream distance can be as- 
similated as an evolution with time t = x/U, using the 
Taylor hypothesis. The spatial resolution of the images 
is 50 pm. Analysis of the scalar spectra revealed no par- 
ticular power law, as well as a cut off at k ~ l^ 1 which 
corresponds to the laser sheet thickness I ~ lrram (20 pix- 
els). It implies that the scalar field measured is filtered 
at this scale I, which is also found to be a typical value 
of the scalar field correlation length. 

c. Mixing from a line source Concentration PDF 
measured at increasing distance from the injector are 
shown in figure[2l It has been compared to the prediction 
of the self-convolution model, shown as a solid line. For 
that purpose we start from the measured "initial" PDF 
at x = 36 cm, and apply solution (f5]), with f(t) fitted 
to the measured variance C2(t). Strikingly, there is good 
agreement with the model based on a pure convolution 
mechanism. We see no tendency for the PDF to approach 
gamma-PDF (dashed line). 
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FIG. 2: Line injector. The experimental PDF at x = 36cm and 
variances at x = 42cm and x = 77cm are used for the prediction 
of the self-convolution model. The fit with Gamma-distribution 
requires only the variance value. 
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FIG. 3: Line injector. Experimental PDF at x = 36cm and 
x = 42cm. The coarse graining is applied to the scalar field at 
x = 36cm. 



In figure [H the initial PDF, pi(xo) has been obtained 
by analysing scalar fields at the resolution of the laser 
sheet. What would happen if the scalar field were coarse- 
grained with a filtering function of length I' > I ? The 
PDF pi'(xo) is expected to become sharper and sharper 
around its mean value, as if the scalar field were "mixed" . 
This filtering process therefore has the qualitative effect 
of a diffusivity. The initial PDF pi(xo) is reproduced in 
figure [3] Let 02(0) be its variance. We then compare the 
coarse-grained PDF pi (xo) with the PDF pi(x) measured 
at a distance x such that their variances have the common 
value C2 = c 2 (0)/2. We also plot the self-convolution of 
the initial PDF pi(xo). The PDF shown are zoomed close 
to the mean value because the filtering process reduces 
the number of independent events for the construction of 
the PDF pi{xq). As expected, the coarse-grained evolu- 
tion is qualitatively close to the temporal evolution, with 
discrepancy for low scalar value. In addition, the self- 
convolution of the initial PDF is in very good agreement 
with the coarse-grained PDF. 

The comparison between the decay of different cumu- 
lants of the experimental PDF pi with the prediction of 
the model reveals a limit of the self-convolution approach. 
On figure 3]), we fit the variance decay with a power law 
C2(x) ~ x~ 2 ' 5 , which is consistent with experimental re- 
sults of The self-convolution model predicts then 
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FIG. 4: Line injector. Evolution of the first cumulants. The 
dashed and dot-dashed line are prediction of the self-convolution 
model, by supposing / ~ x a , where a = 2.5 is obtained from the 
experimental variance decay 
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FIG. 5: Point injector. The fit with the model of pure self- 
convolution is done by taking the first measured PDF (dark grey) 
as the initial condition. 



a decay C3 ~ x~ 5 and C4 ~ x~ 7 ' 5 . We notice fluctua- 
tions of experimental data around this prediction, and 
a slight overestimation of the actual decay of the cumu- 
lants, growing at increasing distance from the injector. 

d. Mixing from a point source The sequence of ex- 
perimental PDF presented in figureOhas a characteristic 
shape that indicates that almost no mixing has occurred. 
There are two sharp peaks : one is the backround flow, 
with a small mixed scalar concentration, and the other 
one corresponds to the existence of unmixed blobs of 
scalar. This second peak persits in the early evolution 
of the distribution. Then the peak disappears, but the 
PDF remains very skewed, showing no convergence to- 
ward a gamma-PDF. It is also clear that the model ((T|) 
does not provide good prediction in this case: there is a 
spurious peak. 

e. Discussion We observe that the self convolution 
model provides good predictions in the case of the line 
injector, but discrepancy clearly appears with the point 
injector. We here discuss the difference between these 
two cases in relation with our initial hypothesis. Since 
the capability of the model to describe the time evolution 
of the PDF depends on the scalar injection in the same 
turbulent flow, the observed descrepancy should not be 
sought in the fluctuations of f(t) (first model hypothesis). 

To discuss the validity of the second hypothesis, i.e. 
the independence of scalar probability functions, let us 



4 



a) 




b) 



poin 


; 




1 


J 


1,1 



512 1024 

x (pixels) 



512 1024 
x (pixels) 



FIG. 6: Typical scalar profile on an image taken at 40 cm from 
the injector. The total length is 5 cm (1024 pixels). The filtering 
length scale 1 (laser sheet thikness) is around 20 pixels. 



consider a one dimensional idealization of the injection 
system: the scalar field is supposed to be a succession of 
voids and colored segments, uncorrelated to each other. 
While the corresponding PDF is a double Dirac function 
in both cases, the PDF of the coarse grained fields may 
be very different. Two coarse grained scalar concentra- 
tions, separated by a distance I, can be considered inde- 
pendent only if I ^> max{Z (i), ^o-(i)}, where .lo(t) and 
l a (t) are typical lengths of voids and colored segments. 
For both the line and the point case, we are in a dilute 
limit: lo(t) ^> l a (t). As seen in figure 6-b, in the case of a 
point injector, the voids length scale 7o is larger than the 
probe length scale I (20 pixels). It has been estimated 
more quantitatively to be around 80 pixels by comput- 
ing the statistics of voids. It is thus not surprising that 
the self-convolution model does not work in that case. 
By contrast, as seen in figure 6-a, it is not possible to 
find void events in the coarse grained scalar field. It in- 
dicates that the typical void length scale is smaller than 
the probe length scale I. 

More generally, those length scales are related to the 
channel and the injector properties. The ratio of their 
cross section surface, respectively ~ rl s /l 2 s and ~ r 2 /I 2 
for the line and the point case, is equal to the mean 
scalar value (a) = l a (t)/ (Jo(t) + I CT (i)). The distance 
r ~ 2 mm is the initial scalar sheet thickness (either the 
width of the vertical tube or the diameter of the small 
cylinder) . This length is then reduced by straining effects 
l a {t) — r I f(t), and the random succession of scalar sheets 
create void segments, while conserving (a). This gives 
an estimation of the void length scale: Jo(t) = l s /f(t) 



for the line injector, and lo(t) = l 2 /(rf(t)) for the point 
injector. As far as r <C l s , we see that void length scales 
differ by an order of magnitude. This is consistent with 
experimental observation: voids events are of the order of 
a few filtering scales I for the point injector, and smaller 
than I for the vertical injector. Note that at sufficiently 
large time, the condition for independence of scalar sheets 
will always be realized. 

/. Conclusion If a gamma PDF is initially created, 
our self-convolution model predicts that it will evolve to- 
ward a Gaussian through a succession of gamma-PDF, 
in exact agreement with Our experiments however 
provide no evidence for a general convergence of scalar 
PDF toward gamma pdf, in contrast with the model pre- 
diction of (llj . Further comparisons of different injection 
systems would be needed to confirm this point. 

In the case of a line injector, our experiments pro- 
vide good support of the self-convolution model. In this 
model, the evolution of the scalar PDF by the straining 
effect is equivalent to coarse-graining. We have directly 
checked this equivalence in the experiments. For practical 
purposes, this self-convolution model has two interesting 
properties: i) it conserves scalar bounds and mean scalar 
value ii) it evolves toward a sharp gaussian at large time, 
by contrast with the LMSE (or IEM) model. 

The validity of the self-convolution model depends on 
the ratio between the cross sections of the injector and 
the channel. This ratio prescribes length scales for the 
fine grained scalar structures (length of colored and void 
segments ), which decrease with time by straining effects. 
The model is correct when the probe length scale I is 
larger than the typical scalar field length scales, to assure 
independence between scalar sheets at distance I. 

The model also relies on the hypothesis of weak fluc- 
tuations for the straining rate. This is typically expected 
for the Batchclor regime of scalar cascade or for moderate 
Reynolds numbers, as in our experiments. Discrepancies 
are expected to arise from intermittency in a well devel- 
opped Kolmogoroff cascade. 
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